Two-stage spin-flop transitions in S = 1/2 antiferromagnetic spin chain BaCu 2 Si207 
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Two-stage spin-flop transitions are observed the in quasi-one-dimensional antiferromagnet, 
BaCu2Si207. A magnetic field applied along the easy axis induces a spin-flop transition at 2.0 T 
followed by a second transition at 4.9 T. The magnetic susceptibility indicates the presence of 
Dzyaloshinskii-Moriya (DM) antisymmetric interactions between the intrachain neighboring spins. 
We discuss a possible mechanism whereby the geometrical competition between DM and interchain 
interactions, as discussed for the two-dimensional antiferromagnet La2CuC>4, causes the two-stage 
spin-flop transitions. 

75.10.Jm, 75.25,+z, 75.40.Cx, 75.50.Ee 



I 

CO 



-a 
c 
o 
o 



> 

00 

m 

o 



C3 



I 

C 

o 
o 



X 



The magnetic long-range order (LRO) of a low- 
dimensional S = 1/2 antiferromagnet (AF) is qualita- 
tively different from that observed in a conventional 
three-dimensional (3D) AF. The absence of LRO in a 
purely one-dimensional (ID) S = 1/2 Heisenberg AF at 
any temperature and a two-dimensional (2D) one at T > 
K was proved several decades ago |Q. Weak but fi- 
nite interactions between chains (ID) or planes (2D) are 
thus indispensable for achieving the magnetically ordered 
state at a finite temperature 0] . In such weakly coupled 
chain or plane systems, the magnitude of Dzyaloshinskii- 
Moriya (DM) antisymmetric interactions, if they exist, 
can be comparable to that of interchain or interplane 
interactions, and as a result, DM interactions cause sev- 
eral behaviors different from those of conventional 3D 
AF. For example, the field-induced gap in magnetic ex- 
citation has been investigated for ID Cu benzoate 0] 
with respect to the effective staggered field induced by 
DM interactions [Q. In this letter, we report exotic two- 
stage spin-flop transitions in the quasi-lD spin system 
BaCu 2 Si207, which are probably due to the competition 
between DM and interchain interactions. The magnetic 
properties of BaCu 2 Af2C>7 (M — Si and Ge) have been 
studied recently |5|-|8| , and it was found that the DM in- 
teractions play a central role in the weak-ferromagnetic 
order of BaCu2Ge20y. Since BaCu2Si2C>7 is isomorphic 
to BaCu2Ge207, we expect the presence of DM inter- 
actions in BaCu 2 Si207 also. The results are discussed 
in comparison with the complicated spin-flop nature of 
La2Cu04, a parent compound of high-T c cuprates ||, 
where the symmetric interplane interactions and the anti- 
symmetric intraplane DM interactions give rise to a com- 
plex spin rotation under magnetic fields. 

A single-crystal sample was prepared by a floating- 
zone (FZ) method as previously reported M . The crystal 
was cut into a rectangle with dimensions of 1x2x1 mm 3 
(axbxc). The orientations of all the surfaces were con- 
firmed by x-ray Laue backscattering. Since the diffrac- 
tion patterns from the ab and be planes are difficult to dis- 



tinguish, we checked the orientation also using an io-W x- 
ray diffractometer. Temperature and magnetic- field de- 
pendences of the magnetization (Mi(T,H), i = a, b, and 
c) and susceptibility (xi(T, H) — ' W — -) were measured 
using a commercial SQUID magnetometer (MPMS-7) up 
to H = 7 T, with which we can control the temperature 
and magnetic field to within ±0.01 K and ±1 Oe of the 
set values, respectively. 

First we show the previously reported low-field sus- 
ceptibility data (H = 0.1 T) in Fig. 0(a) [§. The 
high-temperature behavior obeys the theoretical Bonner- 
Fisher curve [ ]10| |, indicating paramagnetic behavior of 
the uniform ID S — 1/2 Heisenberg AF chain. The 
AF transition occurs at Tn = 9.2 K, below which the 
susceptibilities show anisotropic behavior. % C (T, 0.1 T) 
shows a large decrease below Tn, while Xa(T, 0.1 T) and 
Xb(T, 0.1 T) roughly keep the values at Tn down to 2 K. 
This indicates that the c axis is the principal easy axis, 
as confirmed by neutron diffraction ||g||. The large 
residual susceptibility % c (2.0 K,0.1 T) should also be 
noted. The absence of divergent behavior toward T — 
K indicates that this residual Xc is not the response 
of free spins due to impurities and/or crystal imperfec- 
tions, but rather is intrinsic to the spin susceptibility 
of the BaCu 2 Si207 system. Moreover, the anisotropy is 
also observed above Tn- Two characteristic features are 
found in this region. On the one hand, Xa is larger than 
Xb and Xc at room temperature, which can be attributed 
to the difference in g values roughly determined from the 
oxygen coordination around Cu 2+ ions. In the vicinity of 
Tn, on the other hand, the deviation from the Bonner- 
Fisher curve becomes notable only in Xb and Xc, while Xa 
retains the Bonner-Fisher- like behavior even close to Tn- 
This upward deviation suggests the appearance of an ad- 
ditional contribution to the magnetization only along the 
be plane. 

What characterizes this compound is a peculiar field 
dependence of Xc(T,H) below Tn- As is shown in 
Fig. 0(b), Xc traces three different curves at H = 0.1, 3.0, 



and 6.0 T. Since the 0.1 T field is well below the first spin- 
flop field, as will be mentioned later, \c(T, 0.1 T) repre- 
sents the zero-field susceptibility. Xc(T, 6.0 T) shows be- 
havior similar to that of conventional AF in the spin-flop 
phase, where Xc(T) maintains almost the same value as 
that at Tm- In contrast, Xc(T, 3.0 T) is notable in that 
it follows neither x c (T,0.1 T) nor X c(T, 6.0 T). It de- 
creases smoothly toward T = K, but not as fast as 
X C (T,0.1 T). Tm is almost unchanged and no indication 
of another transition is found. Thus we conclude that 
the three sets of data are intrinsic to BaCu 2 Si 2 07. 

The existence of the three phases is more clearly seen in 
the field dependence of M C (T, H), as shown in Fig. ||(a). 
M c (9.2 K(= Tn),H) has an almost linear relationship 
with the magnetic field up to H — 7.0 T. The data 
at 8.5 K, which is slightly lower than at TV, shows a 
steep magnetization jump at H = 2.0 T, and broad in- 
crease from H = 5.0 T to 5.3 T. Both transitions be- 
come clearer in M c (5.0 K, H). The magnetization shows 
a discontinuous jump at H c \ = 2.0 T and i? c2 = 4.9 T. 
Consequently, it becomes evident that the x c 's shown in 
Fig. |](b) represent the susceptibilities of these three re- 
gions. The phase diagram obtained from the magnetiza- 
tion and the susceptibility is shown in Fig. ||(b) . The low- 
temperature ordered state is separated into three phases; 
for convenience, we call the lower-, middle-, and higher- 
field phases AF, SF1, and SF2, respectively. The phase 
boundaries between AF and SF1 and between SF1 and 
SF2 are almost temperature independent. We also con- 
firmed that the spin-flop transition is absent when the 
field is applied along the a- and 6-axis directions (only 
Xb(5.0 K, H) is shown in Fig. 0(a)). Thus, we can elimi- 
nate the possibility of the mixture of misaligned domains 
being responsible for one of the two transitions. 

The presence of three phases is not expected in a con- 
ventional AF. To construct a proper model in order to 
explain these two-stage spin-flop transitions, we consider 
the following three results. First, the c axis is the prin- 
cipal easy axis at H = T, as was determined by neu- 
tron diffraction |5|,|8j. Second, Xb and Xc are enhanced 
in the vicinity of Tjy. Such behavior is characteristic of 
the AF transition with spin canting, as was observed in 
La 2 Cu0 4 |{n| and Bi 2 Sr 2 Co0 6 . 25 @. In the present 
case, the canted moments are confined to the be plane 
and probably originate from the DM interactions. Third, 
Xc(T, 6.0 T) is almost unchanged below T/y. Such a weak 
temperature dependence is usually observed in a conven- 
tional spin-flop phase, where spins are almost perpen- 
dicular to the applied field with slight canting toward 
the field direction. In a conventional AF without strong 
anisotropy, however, the susceptibility above T/v exhibits 
a normal paramagnetic behavior (Curie- Weiss law for 3D 
and Bonner-Fisher curve for ID), so that the particular 
enhancement of susceptibility is absent. Thus the pres- 
ence of additional susceptibility at 6.0 T below Tjq im- 
plies that spins are more canted toward the c direction 



than expected for only an the external field. The most 
probable source of this canting is DM interactions. 

DM interactions often occur in a spin system with low 
crystallographic symmetry. Dzyaloshinskii first pointed 
out that an asymmetric superexchange interactions such 
as D • Si x Sj is allowed in several antiferromagnets [13[ . 
Later Moriya gave a microscopic basis for the interactions 
Ji"4| |. For the particular case of BaCu 2 Si 2 0y, the intra- 
chain Cu-O-Cu bond actually has a local symmetry (no 
inversion, but a mirror plane including two Cu 2+ and one 
O 2 " ion) that allows DM interactions at this bond. Note 
that DM interactions are perturbative to the main su- 
perexchange interaction. Thus it is sufficient to consider 
DM interactions only between the intrachain neighboring 
Cu 2+ ions. 

Let us now construct a model structure, taking DM 
interactions into account. By considering the local sym- 
metry of the intrachain Cu-O-Cu bond, we can esti- 
mate Dzyaloshinskii vectors as real-space vectors D^ = 
£>((-l) 4 0.86, (-1) 4 0.51,0.07). The magnitude of the c 
component is one order lower than those of the a and b 
components, so that it can be safely ignored. We also 
assume the presence of an easy-axis anisotropy along the 
c-axis as is suggested by Xc(T, 0.1 T). Then the spin ar- 
rangement at H = T is obtained as shown in Fig. 0(a) 
under the experimentally determined inter-chain inter- 
actions: J = 24.1 meV (AF), J a = -0.460 meV (ferro- 
magnetic), J b = 0.200 meV (AF), and 2J 1W = 0.152 meV 
(AF) H. Spin canting and accompanying weak ferromag- 
netic moment per chain are present, but such weak fer- 
romagnetic moments are mutually compensated between 
the a-axis neighboring chains. Therefore, this model is 
consistent with the observation of no spontaneous mag- 
netization. Whenever we discuss the effect of DM in- 
teractions, we should also pay attention to the effect of 
additional pseudo-dipole interactions accompanying DM 
interactions p5|| . However, the effect of this term (so- 
called KSEA interactions) is to erase the anisotropy, and 
under the presence of easy-axis anisotropy along the c 
axis, it does not change the spin configurations shown in 
Fig. ||(a) . The effect of this term will also be discussed 
later. 

In a conventional spin-flop theory, only symmet- 
ric spin-spin interactions are taken into account, and 
therefore, the competition between Zeeman energy and 
anisotropy energy determines the spin-flop field H c . On 
the other hand, DM interactions are antisymmetric, so 
that they can compete with symmetric interchain inter- 
actions under a certain geometry of the spin arrange- 
ment. This idea was first proposed by Thio et al. to ex- 
plain the anomalous spin-flop behavior of La 2 Cu04 [91. 
They discussed the magnetic behavior of La 2 Cu04 with 
respect to a staggered moment and a ferromagnetic mo- 
ment defined as a difference and a sum of two neighboring 
spins belonging to different sublattices. When DM in- 
teractions are sufficiently strong, the spin-flop transition 



occurs in such a way that the weak-ferromagnetic mo- 
ments (initially perpendicular to the field) rotate toward 
the field direction, and as a result, one does not observe 
the discontinuous magnetization jump at the critical field 
that is observed in a conventional AF. The rotation di- 
rection (clockwise or counterclockwise) depends on the 
direction of D vectors. Thus, if the initial direction of 
the weak ferromagnetic moments alternate with every 
layer of La2Cu04, the rotation directions also alternate 
M. Within this rotation process, a competition between 
DM and interchain interactions can occur. 

To demonstrate such a competition, we further sim- 
plify our spin system as follows. First, we ignore the 
effects of Jfc and Jno and take only J a into account, be- 
cause J a is the largest among the interchain interactions. 
Next, we consider the a-axis component of D vectors 
(D a ) only, because D a must be responsible for the ob- 
served oc-plane anisotropy. Then we can focus our dis- 
cussion on four particular spins (black arrows) shown in 
Fig. ||(b), in which white arrows represent D a . Note that 
these spin chains are coupled with ferromagnetic J a . The 
signs of D a alternate along the a axis, which means that 
the initial directions of the weak ferromagnetic moments 
also alternate along the a-axis. When we consider a single 
chain with alternating D a , as shown in Fig. 0(c), mag- 
netic field applied along the c axis causes one spin-flop 
transition at H c . Below the spin-flop field, the weak fer- 
romagnetic moment rotates with increasing field up to 
H Cl which corresponds to each spin also rotating in the 
same direction. 

Now let us take the effect of a-axis neighboring chains 
into account. If the effect of J a is far smaller than that 
of DM interactions, spins may behave in the same way 
as in the case of a single chain. Figure 0(d) shows a 
schematic of the behavior; again we will obtain only one 
spin-flop field. If the magnitude of J a becomes compara- 
ble to that of D a , however, the situation differs. Because 
of the alternating D a , the rotation direction of the weak 
ferromagnetic moment (and thus each spin) is opposite 
to that of the a-axis neighboring chains. In Fig. ||(d), 
a-axis neighboring spins are almost antiparallcl at H = 
H c ; this is unfavorable for large J a . Such an effect seems 
to be negligible at H k, T and H ^> H c , because in 
both situations, spins are approximately parallel to each 
other along the a-axis. Therefore, a new phase can ap- 
pear only around H c . The presence of the SF1 phase 
suggests that the effect of J a is comparable to that of 
D a in BaCu2Si207. In this phase, the a-axis neighboring 
spins must be parallel to each other. There are two pos- 
sible configurations according to the crystal symmetry: 
one is that the spins are pointed approximately in the a 
direction, and the other is that they are pointed approxi- 
mately in the b direction. However, the latter seems to be 
unfavorable, considering the effect of KSEA interactions 
Jl5[ which effectively add an easy-axis anisotropy along 
the direction parallel to D a . This means that spins tend 



to point in the a direction, and thus the former case is 
favored. As a result, we propose a spin configuration of 
the intermediate SF1 phase, as schematically illustrated 
in Fig. 0(e). By increasing magnetic field, spins start 
to rotate within the be plane (AF) up to H c %, switch 
to the ac plane (SF1), and then return to the be plane 
(SF2) at H C 2- This concept is consistent with the ap- 
pearance of additional spin canting in the SF2 phase; 
the alternating D a can produce additional spin canting 
toward the field direction in every chain. Such complex 
spin-flop transitions cannot be observed in a spin sys- 
tem with only symmetric spin-spin interactions, and we 
believe that they are possible only when the symmetric 
and antisymmetric interactions coexist with comparable 
magnitude. 

For further quantitative discussion, we should consider 
the effects of Jf,, Jno, and the b component of D. We in- 
tentionally ignore the possibility that spins continuously 
change their directions from the ac to the be plane. This 
corresponds to assuming the existence of some unknown 
potential barrier between these two configurations. How- 
ever, a recent neutron diffraction study of BaCu2Si20y 
in a magnetic field revealed a sudden change of the spin 
directions, which is consistent with the proposed model 
|l6| . The detailed quantitative analysis will be reported 
in the future. 

In summary, we have observed a new type of spin- 
flop transition, where the low dimensionality and accom- 
panying antisymmetric interaction cooperatively induce 
characteristic spin-flop transitions. The study of DM in- 
teractions is still one of the hot topics of low-dimensional 
quantum antiferromagnets, particularly from a theoret- 
ical viewpoint U%, and we believe that this compound 
will be a well-defined basic material with which to study 
the weakly coupled S = 1/2 spin chain system with the 
DM interaction. 

We appreciate A. Zheludev and E. Ressouche for fruit- 
ful discussions and for sharing their unpublished data of 
spin structures determined by neutron diffraction. We 
also thank T. A. Kaplan, N. Furukawa and K. Kubo for 
valuable discussions and R. Hiwatari and M. Kiuchi for 
technical advice. Work at the University of Tokyo is sup- 
ported in part by a Grant-in-Aid for COE Research by 
the Ministry of Education, Culture, Sports, Science and 
Technology of Japan. 



[1] H. A. Bethe, Z. Phys. 71, 205 (1931); N. D. Mermin and 
H. Wagner, Phys. Rev. Lett. 17, 1133 (1966). 

[2] M. T. Hutchings, E. J. Asmuelsen, G. Shirane, and K. 
Hirakawa, Phys. Rev. 188, 919 (1969); T. Ami, M. K. 
Crawford, R. L. Harlow, Z. R. Wang, D. C. Johnston, Q. 



[3 

K 

[5 

[-: 

[8 

[9; 

[10 

[11 

[12 

[13 

[14 
[15. 



[16 

[17. 



Huang, and R. W. Erwin, Phys. Rev. B 51, 5994 (1995). 

D. C. Dender, P. R. Hammer, D. H. Reich, C. Broholm, 
and G. Aeppli, Phys. Rev. Lett. 79, 1750 (1997). 

M. Oshikawa and I. Affleck, Phys. Rev. Lett. 79, 2883 

(1997). 

I. Tsukada, Y. Sasago, K. Uchinokura, A. Zheludev, S. 

Maslov, G. Shirane, K. Kakurai, and E. Ressouche, Phys. 

Rev. B 60, 6601 (1999). 

A. Zheludev, M. Kenzelman, S. Raymond, E. Ressouche, 

K. Kakurai, T. Masuda, I. Tsukada, K. Uchinokur, and 

A. Wildes, Phys. Rev. Lett. 85, 4799 (2000). 

I. Tsukada, J. Takeya, T. Masuda, and K. Uchinokura, 

Phys. Rev. B 62, R6061 (2000). 

M. Kenzelmann, A. Zheludev, S. Raymond, E. 

Ressouche, T. Masuda, P. Boni, K. Kakurai, I. Tsukada, 

K. Uchinokura, and R. Coldea, Phys. Rev. B (in press), 

also |cond-mat/0012452J . 

T. Thio, C. Y. Chen, B. S. Freer, D. R. Gabbe, H. P. 

Jenssen, M. A. Kastner, P. J. Picone, N. W. Preyer, and 

R. J. Birgeneau, Phys. Rev. B 41, 231 (1990). 

J. C. Bonner and M. E. Fisher, Phys. Rev. 135, A640 

(1964). 

K. Uchinokura, T. Ino, I. Terasaki, and I. Tsukada, Phys- 

ica B 205, 234 (1995). 

J. M. Tarascon, P. F. Miceli, P. Barboux, D. M. Hwang, 

G. W. Hull, M. Giroud, L. H. Greene, Y. LePage, W. 

R. McKinnon, E. Tselepis, G. Pleizier, M. Eibschutz, D. 

A. Neumann, and J. J. Rhyne, Phys. Rev. B 39, 11587 

(1989). 

I. E. Dzyaloshinskii, Zh. Eksp. Teor. Fiz. 32, 1547 (1957) 

[Sov. Phys. JETP 5, 1259 (1957)]. 

T. Moriya, Phys. Rev. 120, 91 (1960). 

T. A. Kaplan, Z. Phys. B 49, 313 (1983); L. Shekhtman, 

O. Entin-Wohlman, and A. Aharony, Phys. Rev. Lett. 

69, 836 (1992); A. Zheludev, S. Maslov, I. Tsukada, I. 

Zaliznyak, L. P. Regnault, T. Masuda, K. Uchinokura, 

R. Erwin, and G. Shirane, Phys. Rev. Lett. 81, 5410 

(1998). 

E. Ressouche and A. Zheludev, unpublished data. 

For example, M. Daniel and R. Amuda, Phys. Rev. B 53, 
R2930 (1996); O. Derzhko and J. Richter, Phys. Rev. B 
59, 100 (1999). 



0.4 



0.3 



0.2 



0.1 



0.0 



1 1 


1 I J-, 


^n 


1 ' 1 


1 1 ' 


(a) f ~ 


: :»- 


" 


E 


,-s* - 


- 


® 25 

b 


- J/P - 

T -l - 


A 


5 20 






5.0 5.5 




D) '^ 


- ;:>j H(i)#- 


" o 8 




J] 


o 






i_ A '° i i i I . . . 


A* / 


"5 1.5 2.0 2.5 y^V I 


h(t) Jy 




Jy 


M 9.2K(H||c) 

J^y • 8.5K(H||c) 
. /X ° 5.0K(H||c). 


-Jy - 5.0 K (H || b) . 

r 



12 3 4 5 6 
Magnetic Field (T) 




6 8 10 

Temperature (K) 



12 




50 1 00 1 50 

Temperature (K) 



5 10 15 20 

Temperature (K) 



FIG. 1. (a) Magnetic susceptibilities along the three 
principal axes around the Neel temperature measured at 
H = 0.1 T. The data were taken from Ref. [5]. (b) 
Low-temperature c-axis susceptibilities measured with three 
different fields. TV is almost unchanged with applied field. 



FIG. 2. (a) Field dependence of the magnetization along 
the c axis at T = 5.0, 8.5, and 9.2 K. At T = 5.0 K, spin-flop 
transitions are observed at H = 2.0 and 4.9 T. The 6-axis 
magnetization shows no anomaly, (b) Phase diagram around 
Tjv. Open and filled circles are the boundaries obtained from 
X — T and M — H measurements, respectively. Solid and 
dotted lines are guides for the eye. 
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FIG. 3. (a) Expected spin arrangement at H — T under 
the effect of DM + KSEA interactions, (b) Four spins on the 
a-axis neighboring chains at H — T. (c) Spin configurations 
of a single chain with DM interactions with the field applied 
along the c axis, (d) Spin-flop process expected when J a is 
negligibly small. Spins are confined within the be plane, and a 
spin-flop transition occurs only at H c . (e) Possible two-stage 
spin-flop transitions. Spins are confined to the be plane in the 
AF and SF2 phases, while they turn to the ac plane in the 
SFI phase. 
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